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Abstract 

In this paper we consider scattering due to a metallic cylinder which 
Is In the field of a wire carrying a periodic current. The aim of the paper 
Is with a knowledge of a far field measurement In between the wire and the cylinder 
to obtain Information such as the location and shape of the cylinder. The 
same analysis Is applicable In acoustics In the situation that the cylinder 
Is a soft wall body and the wire Is a line source. The associated direct 
problem In this situation Is an exterior Dlrlchlet problem for the Helmholtz 
equation In two dimensions. We present an Improved low frequency estimate 
for the solution of this problem using Integral equation methods and our cal- 
culations on Inverse scattering are accurate to this estimate. The farfleld 
measurements are related to the solutions of boundary Integral equations in 
the low frequency situation. These solutions can be expressed In terms of 
mapping function which maps the exterior of the unknown curve onto the exterior 
of a unit disk. The coefficients of the Laurent expansion of the conformal 
transformations can be related to the far field coefficients. The first far- 
fleld coefflqlent leads to the calculation of the distance between the source 
and the cylinder. The other coefficients are determined by placing the source 
In a different location and using the corresponding new farfleld measurements. 
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Introduction 


In tills note we are Interested In the scattering due to a inetaUlc 
cylinder which la In a field caused by a wire carrying a periodic current. 
The goal Is to obtain Information on the location and shape of the scatterer 
from a knowledge of the far field. A similar problem but with plane wave 
Incidence has been recf’'ntly studied by Colton [1] and Colton and Klelnniann 
[2]. A rigorous analysis on forward scattering In this situation has been 
recently given by Harlharan and MacCamy [3]. If the metal Is not made of 
ferromagnetic material then one will obtain an Interior-exterior problem 
with the boundary of the conductor as an Interface. However, for a ferro- 
magnetic conductor one will obtain an exterior Dirlchlet problem. The de- 
tails of this discussion are available In [3]. The analysis In what follows 
Is also applicable to acoustic wave propagation In the situation that the 
scatterer Is a soft material. The precise model of Interest Is Illustrated 
In Figures 1 and 2. 
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Let f2 be the cross section of the cylinder in the x-y plane, Q 
be Its exterior and F be Its boundary. The wire which Is located at 


carries a periodic current of the form Re(I 



Figure 2 

The derivation of the associated general problem is discussed in [3]. For 
details we refer the reader to this paper and present here the particular 
problem. Find u € U C^(F) such that 

! u - Up e C^(fi+) U G°(D, 

A u + 3^u * 0 in 

u = 0 on r , 

u - Up is the Sommerfeld's radiation condition. 

Here u„ is the potential due to the wire and is prescribed by 

p 

“S " - f «0^’^3l2L-2?ol> (1-1) 

and 3 is a non-dimensional parameter which depends on the frequency (a. 

For moderate frequencies such as household 'urrents or AM transmissions the 

—9 -2 

value of B ranges from 10 < 3 < 10 . This is the size of 3 which is 

2 

of interest in this note. Since 3 is extremely small one might take 3 “ 0 
in problem 1, but there are two difficulties. First, the radiation condition 
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will no longer be valid and second the Incident field potential (1.1) 
will become logarithmically Infinite, However, It turns out that the 
following theorem Is true. 


Theorem 1.1 


u ■ V + 0 (1/log 3) 

where v satisfies the following problem 


as 3 0 


II 


where 



V 6 c°(r) 

- Vq € C^(fi'^) U C°(D, 

Av ■ 0 In 

V *• 0 on r , 

bounded as |)c| «> 


(x € m^) , 


log|x-XQ 


(x€ . 


Thus the solution of problem II will be the limiting value of thfi solution 
of problem I as 3 0. 

In order to solve the Inverse scattering problem several details of 
forward scattering are required. In this process we shall prove theorem 
1.1 and also provide an Improved estimate: 

Theorem 1.2 

CF 

“ ■ W- 

where C C IR , € C- and fQ(x) € C^(f) Is a uniquely determined function. 

Hence our results on inverse scattering would allow one to obtain the 

2 

shape of the boundary for low frequencies accurate to 0(3 log 3). 


-4 


2. Forward Scattaring 

Hera we shall Indlcace the solution procedure for problems 1 and II 
If the boundary, F (sufficiently smooth), is given. Wa shall use simple 
layer potentials for the construction of our solutions. 

Let 

/^<y(z)G5<x.x)d8 » (<S « m ) , 

r 

log|x- y| If « - 0, 

(-f if Mo. 

Than the following properties hold: 

(1) ^ 0 in 3R^ 

(11) € c(m^) 

(111) satisfies Somraerfalds radiation condition If 6 0. 

(Iv) l/g a* (^ff(x)da^)log|x| as |x| If «S « 0. 

Notation: 

[cf] (x) - F^[o](x) for X € r. 

We shall use this notation throughout this section. 

Thus one can seek a solution of problem I in the form 

* *^gt<j>Kx) + »0(x) (x * ^5^)» (2.1) 

Then (2.1) will be a solution of (I) provided 

III . + Uj^ « 0 on r. 

Similarly, we seek the solution of problem (II) in the form of a simple layer 
except incorporating a constant C € 3R which is: 


where 

G^(x,x) - 


( 2 . 2 ) 


v(x) " ^ (X ® 


Demanding that v satisfies the boundary condition and the condition at 
infinity, we have, 


IV 


! 


+ Vq + C ■ 0 


C ijj(s)ds " -It 

•^r 


on r 


Remark 2.i 

(a) The constant C is incorporated in (2n2) in order to make the system 
of integral equations (IV) consistent. This is a standard argument. 

(b) The uniqueness of problems X and IX are also standard as is the exis- 
tence of <J) and (>|^C), see for example [3], [4). Our low frequency theory 
is based on the exiBtence proofs presented in fS]. 

To present our main result we need several lemmas. We shall begin by 
stating them. (The proofs may be found in [3J.) 


Lemma 2.1 

There exists a unique solution (f^jF^) C CfF) x IR for the system of 
integral equations 

- ^0 

^ (2.3) 

JJ^fQ(s)d8 - 1 . 

Moreover, the solution is given by 

^0 " '■■^lr^°8|f(z)|, To log|f(z)| (z C D (2.4) 

where w - f(z) Is the unique mapping function which maps F on to the ex- 
terior of a unit disk ju] « R such that f(z) z as z 


" 6 " 


Renmrk 2.2 

If we nmp P onto the exterior of a unit disk with f(z) 'V az as 

z « (a > 0, a”^ is the mapping radius) then |w| ■ |f(z)I • 1. Thus 

Fo-O. 

As an aside we would like to point out that the result of this lemma 

enabled us in [3] to prove the existence of the system without the condi- 

tion that the mapping radius should not equal one. This is different than 
the proof available in [4]. The following two lemmas may also be found in 
[3]. A version of these proofs with less regularity assumptions are avail- 
able in [5]. 


Lemma 2.2 


There exists a kernel P(s,g) of the form 


P(s,0) " - ~ cot + R(s,a) , 


with R(*»*) analytic such that 

f(s) “ " J P(s,a)h' (a)dg , 

satisfies 

%EfHs) - h(s) -f r(h), 

where F(h) (a constant) is a functional of h given by 

r(h) 




2tt 

Remarks 2.3 

Let us observe for later purposes that when h = 1 


(2.5) 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


and 


r(i) - -1, 


(2.9) 


fp - 0. 


-7- 


It Is shown In [3] the coinblnatlibn of Lemma 2.1 and 2.2 yields 
the existence of (III). In particular when we set h(s) - -Vq(s) we 
have the following 


Lemma 2.3 

The unique solution of the system IV Is given by 


ij) - fp + ctfQ, (2.10) 

C - - r(h) - afQ, (2.11) 

where a € » Is a functional of h given by 


a(h) - “ 1 “ / fp(s)d8. 


( 2 . 12 ) 


Remark 2.4 

In Lemma 2.2 the normalization of T over [0,2 it] is used for 
convenience. 

Now we return to integral equation III. First of all by the expansion 
of the Hankel function (z) for small z we have 

- i H^^^(z) « jlog(z/2)+ (y- (TT/2)l)j + O(z^logz) . 


Using this equation (III) for small 3 leads to the following equation; 


t^o[(<>] + >" / <}>(8)ds + Vq + m + 0(3^1og 3) - 0, (2.13) 

where 

2^ +Y-|l) . 

In (2.13) the term 0(3^1og 3) contains Integral operators involving ({). 
Motivated by (2.13) we consider an associated equation that 


8 - 


m* 

where 


+ m y *f n * h> 


h * -Vq. 


It will be argued later that 

tj) - $ + 0(0^ log 0). 

Leinma 2.4 

The solution of (III*) Is given by 


Cf„ 

^ * + 


Proof 

This Is done In two steps. First we seek a solution of III* In the 

form 


^ " ""'®1 ® 2 ’ 

where g^^ and g 2 satisfy 

i^otSll + “ (2.14) 

^0^®2^ + ™ Jj^g 2 <s)ds - h. (2.15) 

We look for the solution of (2.15) In the form 

82 - fp + XfQ, a^C). (2.16) 

Demanding 82 to satisfy (2.15) we have 

and using lemmas 2.1 and 2.2 we have 

-r(h) “ m fp(s)ds 

X - X(h) . (2.17) 

m + Fq 

Thus 82 Is known. The solution of (2.14) Is obtained in the same way 
setting h 1 . 
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gj^ " fp + X(1 )£q. (2.18) 

We observe by remark (2.3) that 

X(l)«;p^, (2.19) 

«• X(1)£q, 


and hence 


~r(h) - (1+ ^fp(B)ds)m 

I - fp (X(h) - mX(l))fQ - £p + * STTT- ^0 • 


Using the definition of a(h) In (2.12) we have 

-rfh) r^g ^ (m4-rQ)tt 


? - £, . 


m + r. 


fn» 


or 


r(h) 4- r^a 

$ w fp + afQ fjj, 


Now using (2,10) and (2.11) we have 


C£, 




and the lemma. 

Remark 2.5 

Uniqueness of the solution III* can be argued by considering the homo- 
geneous equation of 111* 

Now we shall return to our main theorems 1.1 and 1.2. We shall describe 
the proof of theorem 1.2 which also proves theorem 1.1. First of all we Indi- 
cate the argument that 

({>-$ + o(0^iog a). 


( 2 . 20 ) 


c.« mu:..: t*s 


1 


f]0'* 


To do this wo simply consldor tho dlfforenco oC oquotion (III) and (III*) 
which is 

+ m f (f})-|)(@)ds 0(g^log|) , 

r 

The operator on the loft hand side is invertible • This can be shown in 

the same way we did for equation (2.15) with h ■ 0($^log 6), and the 

inverse operations are performed on terms of order 0(6 log $), Thus we 

have n 

(|) * <|) i- 0(6^og 6). 

Now we return to the simple layer representation (2.1) 

u(x) * ^gt<|)J(2i) + 

From the expansion of Hankel function for small 6 we have 

wCx) “ ^ + m + Vq + 0(6^1og 6)« 

Substituting the relation (2,20) we have 

u(x) * Uq(^J(x) + m^$J(8)ds + m + Vq + 0(6^1og 6), 
and hence from lemma 2.4 we have 

0 j, 

JJfQ(s)ds + m + Vq(x) + 0(6^’log 6)* 


But 


I i|j(s)ds - -1, 

r 


/fQ(s)ds 

r 


1 . 


>5? 






and 


1 


“ 11 " 


Thus, 

u(x) - ■*” ^0^^ ■*■ OC0^1°8 6)» 

Observing r 

Cm _ ^0^ 

^STTq ^ ” fflTl^ » 

we have 

»(*) " ■** ° **■ 0 )- 

From the simple layer potential for problem (II) we have 


U(x) -v(x)+^ jUot£o](x)"ro 


+ 0 ( 0*^108 0 ). 


( 2 . 21 ) 


This proves theorem 1.2 and theorem 1.1 follows immediately. 


3. Far Field Measurements and Inverse Scattering 

In this section we follow the work appearing in [1] and [2J. The main 
idea is to obtain a sequence of moment problems relating the Fourier coeffl" 
dents of the farfleld to the coefficients of the Laurent expansions of the 
conformal transformation which maps the exterior of a circle onto the exterior 
of t^e unknown curve. As in [IJ and [2] it turns out that these moment pro- 
blems are solvable in terms of the mappivig radius a, which in turn is deter- 
mined by placing the source In a different position. The main difference here 
Is that the location of the source can be calculated, thus one Indeed finds 
the location of the scatterer. Such a result is useful in radar detection, 
where one Is Interested more in the location of the radar than the shape. 

The physical situation is portrayed in Figure 3. To avoid repetition 
we shall cite [1] or [2] wherever possible. 
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Figure 3 

We assume that the far field measurements are available on a circle of radius 
r and let the source be at a distance jx^| such that |>^| > r. From pro- 
blem I the scattered part of u is given by 


u®(x) “ - f/<l)(z)H^^^|3|x-xl)ds . 


0 

From the asymptotic behavior of the Hankel function one has 

u®(r,0) 




TTPlr 


F(0;p) + 0 


Cl) 


(3.1) 


where 


FOiS) . - . 


(3.2) 


F is known as the far field. Expanding F(0;p) in a Fourier series we have 

F(0j3) » 2 an(3)e^"®. 


n» -00 


where 
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1 

2tt 


-tr * 


-ln0 - 10P cos 0-(P 


ds^dO 



(3.3) 


Reciark 3.1 

Usually one ineasures a finite number of We shall be concerned 

only with a^ for n ^ 0. Negatively Indexed coefficients do not play a 
role In our calculation (see [2]). 

The Taylor series expansion of Bessel's function In (3.3) yields 

a (3) « -JL^ ^(jr)p"e"^"'^ds . 

" 2"n! •' ^ 

Define 

2”nl a,(6)i" 

u (3) S . (3,4) 

n ^n 

Therefore 

M„(3) - - /5>(y:)p”e"^”‘^d8 . 

r 

Now using (2.20) we have 

- - /?i(y)pV^"‘^’d8 (mod(0(3^ log 3)). 

P JL 

2 

We shall drop mod(0(3 log 3)) for convenience of writing. Observe that 
from Lemma 2.4 we have 

^<x) “ >l^(z) + (3*5) 


and from the simple layer potential (2.2), i|j(y) Is given by 

i>(^) - - |^(v-Vo)(y.), r (3.6) 


where the normal derivative Is In the outward direction. We want to map 
r onto the exterior of a unit circle by 


w ■ f (z) ■ az + 0(1) , 
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wheire a is the mapping iraclius. Thus Iq » Q. We would also like to obtain 
g(z) In terms o£ this mapping function f(K)i For this we have; 

Lemma 3.1 

Let w « f(z) " az + 0(1) bo the unique mapping function which maps 
r onto the exterior of the unit circle |w| •* 1. I'hen the solution of 
problem II is given by 

v(z) log|f(z) - f(zQ>l - lQg|f(z)f (Kq) “ ll » (3.7) 

where Zq is the location of the source and the bar denotes conjugation. 


Proof 

We make use of the uniqueness of the solution of problem II. First 

observe that the expression given by (3.7) 

(i) satisfies Laplace's eqviation. 

(ii) is bounded as z (since f(z) a z as z *^ «). and 

(lii) satisfies v(z) * 0 for z # T (this follows from the fact that 

■ r ' 7 " '- r ■ f(z) for z C r) . 
f (z) 

(iv) furthermore, v-Vq is regular (since log] f (z) - f (Zq) j has leading 
singularity log | z - Zq | ) • 


We now consider the behavior of v(z) - ^0 First 

observe that from the Taylor series of loglz-Zgl for jz^j large that 


Vo(») - 


X 

log|z- ZqI k loglZf,! . 


Since |zq| is large £(zq) » szq + b where b has the form 

h h 

° ^0 "O 
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We want to use this in the exptoaalon (3.7) for v. First of all 
logjfCz) - f (Z q)| - 


>\t 

where R.P denotes the real part. Since a>0, f(zQ) " H>®r®fore 

the second expression in (3.7) becomes 

logjf (z)f (Zf.) - l| « log| (az^ + B)! (z) - l| " loglaz^f (z) j-t log] |l- b 

^ ^ { a5Q\f(z) 

j Logjaz-xj ■4-log|f(z) I - i b\L 

( " Uq f(z) /j 

Tims noting R.p(';;^) * R.P(-|) we have 

“o 

V<«) - i log|£(z) I - ^ (lia.i £(z)) . (J.8) 

0 Zf, 

We now consider two special cases. 

Case (i) Zq 5 (a,0) “ a (Source on x axis) 

In this case 

V(z) - - loglf(z)l -*2J^^^-R.P(f(z)-;^). 

Case (11) Zq (0,a) - id (Source on y axis) 

Then 

v(z) - - Ip log|f(z)| -^ 3 -I.P(f(z)+j^) , 


log|l:(z) - (aZQ+b)l 

j fl S5 ^ I 

0 aziQ 

loglaz^l + log|(l- 
logjazQl - R.p(liS^l) , 


where l.P denotes imaginary part. Therefore, 


(3.6) gives 


to. 

“ l0Bjf(z)l - I.P(f(z)+-|^) for ZQ-ia. 

(3.9) 

that Zq 1 “ ct. In this case $ in (3.5) gives 

|j(f(o)- -rfe) + <S+t)i|jl08|£<0)|. (3.10) 

Wo uoo tills in the oppression for ohloh la 

M “ ~ f $(z)P^e 
n oTp 

Taking complex conjugates 

W *“ - / |(z)p"e^”‘*’<is, 

n «fp 

or using z * pe^*^ and (3.10) we have; 

^n ■ - lib /j, -fb’""''” 

- ^ (|+1) jr|j log|f(z)|z''ds^ (3.11) 

To deal with the integrals in (3.11) let us make the following observation: 
IjR.PCf(r)-j^) .-ljMf(r)-j^), 
by the Cauchy-Rieraann equations. 

®bar denotes conjugation. We used In (3.11) the expression that 
^ since (fi and Iq are real valued functions. 



First we consider the case 


1 
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But 


Thus 


Similarly, 


Since 


and 


R.P(f (z) - » 0 for z € r. 

1^ log|f(z)| -- 1^ arg f(z). 

0) ■ f(z), 
z - f"^(w), 


we can rewrite (3.11) as 


^n“ 2 ^ J)[f“^a))fd8 + ^(|+1) /|^ arg w[f"^(u»)]"d8. 


(3.12) 


Observing that 


9 9 i + 

argo) ds - gjjj argo) dw dw on |o)| - 1, 


we have that (3.12) can be written In the form 

Vn - / I 

|w| - 1 m w 


- 1 


i [f“^(a))]"dw. 


(3.13) 


Since f(z) has the Laurent expansion 


see [2] for these calculations. 


1 
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f (z) " az Bq + ~ + I ... I •• w, € C, i *- 0»l,2|. . . ,) 

tSt 

••1 

we have f has a Laurent expansion of the form 


£'1(U1) - f+ l^o +-“ + -!+ .•■■ 


(b^C J, i - 0,1,2,....) (3.14) 


Since F Is smooth this series is uniformly convergent on |w| - 1, We 
substitute (3.14) in (3.13) and evaluate the residues for each n. We 
shall indicate the first three and the rest can be done recursively. For 
simplicity we shall take Bq » 0, since it represents Just a translation. 


For n * 0 


Hq !• 

U ~ 


(3,15) 


For n ■ 1 


V' 


(3.16) 


For n “ 2 


^“2 


U . 
m ®2 


m 



(3.17) 


etc. From (3.16) one can determine b^^ in terras of a and c»; knowledge of 
bj^ yields bg in terms of a and Ct. In a similar way we can obtain the 
higher order coefficients. The most Important feature is from equation (i) 

C can be determined. 


l.e. , C ■ m(lJQ - 1) . 


(3.18) 


We shall later show that 


C 







Thus (3,15) indeed measures the locations of the source and assigns origin. 


4* 

' This translation will apply for both the origin of the scatterer and the 
location of the source. 
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Now knowing the origin we place the source at Zq - ia , l.e,, on the 
y axis. [See Figure 4]. 



Figure 4 



20 


but 

m 0 jfor z C r. 


Therefore 




(3.ai) 


and when we substitute (3.21) Into (3.20) we have, 




n 


^ / l;(“+s)iriwi"d., + i(§ 4 .« /A(rito)]"d«. 0.2?.) 
nul-1 m •' 


Computing the residues In (3.22) as before we have, 


IlJ - “ (?+l) , 

m 

^ if 1 . . >. 

I'l ■ s'- ’’!> • 


(3.23) 

(3.24) 


Thus (3.16) and (3.24) are two equations for a and In terms of known 
a, and Adding these equations we have 

q2 „ ^ ^3 25) 

a(fll + [lj^) 


Thus the mapping radius a can be determined. 

Remark 3.2 

Since In practice only a finite number of a^'s can be measured we can 
determine only a finite number of b^^'s in (3.14). In [2] an error estimate 
is given when (3.14) la truncated after N terras. This applies in our case 
too. 

Remark 3.3 

In determining a we used the fact that that 


^ “ 27T I ^0 I ' 


If we consider the simple layer potential representation for the solution 
of problem II 

v(x) ■ ^ ^>i'(y)iog|*-i|d8^ ^ iobIx-XqI 

for |x| < |xq|* (x * Fourier expansion ^?f log|x“xl gives no con- 
stant term In the Integral. However 

i log|x-X(|| - If logU- Z(,| - ^ log|zg| + ^ <^)" - 

which has a leading constant log|zQ|. Sldce v(x) “ 0 for jc€ F, 

C + ~ logjzQl must be zero. 

Remark 3.4 

In our Laurent expansion (3.14) we took b^ ■ 0. This was possible 
since the corresponding translation In the z plane applied both to the 
Bcatterer and the source by the same amount. However, one can In fact 
calculate bp from equations (3.13) and (3.20), [6]. 
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